Dynamics of spin transport in voltage-biased Josephson junctions 
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We investigate spin transport in voltage-biased spin-active Josephson junctions. The interplay of 
spin filtering, spin mixing, and multiple Andreev reflection leads to nonlinear voltage dependence 
of the dc and ac spin current. We compute the voltage characteristics of the spin current (/s) 
for superconductor-ferromagnet-superconductor (SFS) Josephson junctions. The sub-harmonic gap 
structure of Is{V) is shown to be sensitive to the degree of spin mixing generated by the ferromagnetic 
interface, and exhibits a pronounced even-odd effect associated with spin transport during multiple 
Andreev reflection processes. For strong spin mixing both the magnitude and the direction of the 
dc spin current can be sensitively controlled by the bias voltage. 

PACS numbers: 72.25.-b, 72.25.Mk, 74.50.-|-r 



By controlling the quasiparticle spin degree of freedom 
in solid-state heterostructures, spintronics provides a new 
paradigm and possibilities for making devices and circuits 
with improved functionality over conventional electron- 
ics Although spintronics research has been mostly 
focused on hybrid structures of normal metal (N), semi- 
conductors, and ferromagnetic materials, superconduc- 
tors has increasingly proved to be a useful component 
in constructing spintronic devices. This is because the 
proximity coupling of superconductor (S) to ferromag- 
net (F), on the scale ranging from nanometers to mi- 
crons, offers a unique opportunity to probe the spin de- 
gree of freedom and coherently control the spin polar- 
ized current 2, 3]. For example, two leading methods 
of measuring the spin polarization of ferromagnetic ma- 
terials, spin-dependent tunnelling f3| and point contact 
Andreev reflection spectroscopy 5], both rely on mea- 
suring the conductance of superconductor-ferromagnet 
junctions [H, 01- Recently spin-triplet supercurrent was 
created to flow through the superconductor-ferromagnet- 
superconductor (SFS) Josephson junctions 0|- The spin 
current transported through in SFNFS structures was 
also predicted to have properties with the potential of 
making superconducting spin-transfer devices 0, . 

The equilibrium properties of SF hybrid systems have 
been extensively studied in recent years 0, [sj. However, 
nonequilibrium spin-transport properties, which are es- 
sential for the purpose of manipulating spin current flow, 
are far from being well understood. As an example, the 
SFS TT junction has been proposed for constructing an 
environmentally decoupled quantum bit [10]. Optimiz- 
ing its performance requires knowledge of its dissipa- 
tive dynamics under nonequilibrium conditions, such as 
the current-voltage characteristics which have only been 
studied in certain limiting cases 



11 , 12, l3| , our method allows one to find analytically the 
time-dependent nonequilibrium (Keldysh) Green's func- 
tions near spin-active interfaces. We apply the formal- 
ism to compute the spin current-voltage characteristics 
of voltage-biased SFS junctions to show that SFS struc- 
tures can serve as highly tunable spin filters as well as 
spin-current oscillators for spintronic applications. 

We consider SFS point contacts to illustrate the dynam- 
ical properties of spin-active Josephson junctions. The 
left (1) and right (2) leads are assumed to be spin- 
singlet, s-wave superconductors (S) in the clean limit. 
The ferromagnetic metallic layer (F) is assumed to have 
a homogeneous magnetization fixed in the z direction by 
anisotropy energy. We also assume the F layer is spec- 
ular and that its thickness is much less than the super- 
conducting coherence length and the mean free path in 
the leads. The contact is characterized by an S-matrix 
[l^ . Each S-matrix component is a 2 x 2 spin matrix, 
which for a fixed quantization axis for spin, z, is diagonal: 
Su = S22 = diag (/Rfe"'/^, /RTe"''^/') , ^21 = ^12 = 
diag (iy^i>i"e"'/^, i^/Oje"*''/^). Da is the transmission 
probability for normal state conduction electrons with 
spin projection a —'\, j, and Ra = 1 — Da- In general 
D^ ^ Di, so the F layer acts like a "spin filter"; i}, re- 
ferred to as the "spin mixing angle" , is the relative phase 
between the transmitted (and reflected) spin up and 
down electrons. Spin mixing at the F layer leads to ro- 
tation of the direction of the spin polarization analo gou s 



In this paper we present a new, general theoretical frame- 
work to compute the time-dependent transport proper- 
ties of SF hybrids. Compared to previous approaches 



to Faraday rotation of linearly polarized photons 14 1. 
Particle-hole symmetry relates the S-matrix for quasi- 
holes to that for quasi-particles, jgij(p/||) = 5ji(— p/n )*'' 
[20| . which simplifies for junctions with inversion sym- 
metry to S_ij — Sij. Da and 1} can be calculated from 
material parameters such as the exchange field, the Fermi 
wave vector mismatch (FWM), and the F layer thickness; 
see section III of Ref. JLOi for detailed discussions on the 
interface S-matrix. We emphasize that FWM plays an 
important role in FS junctions. For example, as shown in 
0, the subgap conductance of a FS junction is very sen- 
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sitive to FWM and a nonmonotonic function of spin po- 
larization. For transition metal ferromagnets such as Co 
and Ni, the transmission probabilities, D^, range from 
0.2 to 0.9, while d can take any value between and tt, 
depending on the F layer thickness [2l[. 

For NFN junctions with normal metallic leads (N) under 
voltage bias V, the spin current flowing through the junc- 
tion is due to spin filtering and linear in V, /s = Gs V. 
The spin conductance Gs oc (D^ — Di) is independent 
of the spin-mixing angle. Spin mixing, on the other 
hand, plays a central role in transport through SFS junc- 
tions. For example, it can give rise to 7r-junction behavior 
[1) [iBl > '^^11 complex features in the IV character- 
istics of the charge current, i.e. the "sub-harmonic gap 
structure" (SGS) For voltage-biased SFS contacts, 
the effective transmission probability for spin up (down) 
quasiparticles in the superconducting leads depends not 
only on D-f and D^, but also on the spin-mixing angle 
I'd) and the time-dependent Josephson phase difference 
across the contact, (j>{t) — 2eVt/h. The latter leads 
to inelastic scattering of quasiparticles of energy e by 
the spin-active interface, i.e. multiple Andreev reflection 
(MAR) [3], which transfers quasiparticles inelasti- 
cally into side-band states with energies e„ = e -I- nhuj, 
with hio = eV. As a result a time-dependent spin current 
flows across the junction which has dc and ac components 
with frequencies 2na;, 

oo 

Is(t) =Io + ^(I^, cos2nwt + I^sin2nwt). (1) 

n=l 

Our goal is to flnd the dependence of {lojln^''} on the 
bias voltage and scattering parameters {Z^q,??}. 

Our method for computing time-dependent transport 
properties is based on formulating the quasiclassical 
equations of nonequilibrium superconductivity [2^ in 
terms of particle-hole coherence functions and distribu- 
tions functions, collectively known as the Riccati ampli- 
tudes. These functions obey Riccati-type transport equa- 
tions and uniquely determine the quasiclassical Green's 
functions. The Riccati formulation is discussed in Refs. 
I9L 23l |: we follow the notation of these authors. For su- 
perconductors in the clean limit the coherence function 
7^ (7^) is the local probability amplitude for a hole (elec- 
tron) being converted into an electron (hole). The distri- 
bution functions and describe the nonequilibrium 
occupation of the particle and hole states. Quasiparticle 
scattering near the contact, including spin filtering, spin 
mixing and MAR, is described by a set of time-dependent 
boundary conditions for the Riccati amplitudes. We have 
derived these boundary conditions for spin-active inter- 
faces in Ref. [13] ■ Here we present exact solutions for 
the boundary conditions for the case of a constant bias 
voltage. The local Green's functions at the contact are 
constructed from the Riccati amplitudes and the spin 




FIG. 1: Diagrammatic representation of Eqs. ((2}. The solid 
(dash) arrows represent electrons (holes) , the vertical line in- 
dicates the interface. Branch conversion amplitudes (7^ and 
7^) are shown as semi-circles connecting electron and hole 
lines, interface scattering amplitudes (Sij) are shown as circles 
(o), and the effective reflection and transmission amplitudes 
{r^i and t^i) are indicated by bullets (•). Only the leading 
order multiple scattering processes for r^i are shown. 

current is calculated from the Keldysh Green's function, 
2^ . We have checked that our method reproduces known 
results for nonequilibrium charge transport through spin- 
independent Josephson junctions |l7i . ilsL l24i. i25|. as well 
as the SFS structure studied in Ref. [Iq . 

We set the phase = and <f>2 = 0(i) and work in a 
gauge where (j){t) = 2eVt/h. The corresponding Riccati 
amplitudes for the two leads are (in the energy domain): 
7f (e',c) = 27r<5(e'-6)7f(e), xf (6',e) = 2^<5(e' - c)xf (e), 
72^(e',e) - 2TTSie' - e + 2u;h^{e - cu), and (e',e) = 
2TTS{e' — e)x2 + t^)- The remaining Riccati amplitudes 
are obtained from those above by symmetry relations: 
7/^ = [7f]^andg(p/,e',e) = g*(-p;,-e',-e) Q. Note 
that ■y^ and 7I*' act as ladder operators in energy space. 

Consider a set of scattering trajectories {li, lo, 2i, 2o} 
near the contact ( "i" =incoming and "o"=outgoing), c.f. 
Fig. 1 of Ref. [19(1 . Boundary conditions relate "outgo- 
ing" and "incoming" Riccati amplitudes via the interface 
S-matrix. For example, on trajectory li is given by. 



rf = r^^o7«^l,-ft« 07^^12, (2) 

where = 8}^ — lfs}ij ° 7,f , and the convolution is de- 
fined as [x o ?/](e', e) = / de"x{e' , e")y{e", e)/2n in energy 
domain. Equations ([2]) are the summation of multiple 
scattering amplitudes represented by the diagrams shown 
in Fig. [T] Constructing from the Ricatti amplitudes 
then leads us to evaluate the product, 

A EE (1-rf o7f)-ior^^ - [/3ii-/?2io/32l,io/3i2]^^ (3) 

The ladder properties of 7^ and 7^ lead to A(e', e) = 
27r(5(e' — e — mu;)Am(e), where Am(e) has the phys- 
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ical interpretation of the transition amplitude for qusi- 
particles of energy e in lead 1 to be scattered into energy 
The calculation of {A„j} requires inverting an infi- 
nite dimensional tri-diagonal matrix in energy space, or 
equivalcntly solving the three-term recurrence equation, 



where the coefficients are defined by 

Om — /3ll(e„j) - S']j/922^(em+l)S'|2 

Pni = 7? (em)S;^i72^(em-l)/322^(em-l)S'|2 , 
<?m = '5'2l/322^(<^m+l)72*''('^m+l)i5i27f (em+2) ■ 



(4) 



(5) 
(6) 
(7) 



Equation ^ is solved by introducing transfer matrices, 
Z^, such that A,„ — Z+A„j_2 for m > and = 
for m < 0. The can be expressed as matrix 
continued fractions Z+ = [om + qmZ:^^_f_2\^^i~Prn), Z:^ = 
[om + PmZZn]~^{—qm), and evaluated using standard 
procedures |26j . The resulting Z^ allow us to compute 
Am for even m starting from Aq — [oo+poZZ2+QoZ2]^^ ■ 
The amplitudes for odd m vanish identically. 



The angle- and energy-resolved spin-current evaluated 
to the left of the contact is given by Ji(p/) 
g{pf)TT [(TT3{q^^ — Qig)] /2ni, where the trace is over 
particle-hole and spin space, g{pf ) = NfA^{vf ■ h), Nf 
is the density of states at the Fermi level, v/ is the Fermi 
velocity, A is the contact area, and n is the normal to the 
junction interface. We evaluate the current by expressing 
T'aiSii —01^) in terms of the local Riccati amplitudes and 
carry out the trace over particle-hole space. We obtain 
a sum of 8 terms; e.g. the contributions from electron- 
like quasiparticles from lead 1, which are proportional to 
, are A o xf o A''' and 7^" o A o o A'f o 75*. Simi- 
larly, contributions from and if are expressed 
as convolutions with transition amplitudes B, C, D, and 
E, which can all be expressed in terms of A [see Eq. ^]. 
We evaluate the convolution in the energy domain, then 
Fourier transform to obtain the trajectory-resolved dy- 
namical spin current, 

Ji(p/;i)=.9(p/)Re^e*2fe.*y' ^Sp [^K„fe(e)] , (8) 

k,m 

where Sp is a trace over spin space. The kernel Kmfe(e) 
consists of 8 terms which describe interference between 
inelastic scattering e — > e2m and e —^ ^2(m+k) '■ 

Kmfe(e) = H™fe(e) -I- 7f (e2m)H„fc(e)[7f (e2{m+fe))]^ 
H™fc(e) = Aa^xf (e)Atj^^^j -l-Bama;^ (ei)Btj^_^^j 

- C2mi2^(e-l)Cj,(^^j.j, 

B„(e) = [Am-2(e2)7f(e2)£L72''(ei)-A„(6)Stj/32-/(6i), 
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FIG. 2: The dc spin current (in units of GsA/e) as a function 
of voltage bias V at zero temperature, for different transmis- 
sion probabilities, D| and Di, and spin-mixing angles, i9 = 
(blue), 7r/2 (green), vr (red). 



C™(e) = {D„(e)S:i2 - [A™+2(e-2)Sti 

+ B,„+2(e-2)5'j2]5i27?(e-i)}/322^(e-i), 

D„(e) = A„(e)7f(e)Sli+B„+2(e-2)72''(e-i)5l2, 

E„(e) = U^„,)[A„,ie)Sl+B^ie)Sl^]jZ\e). (9) 



Thus, the calculation of the spin current is reduced to 
the calculation of the {A„i}. Integration over all con- 
tributing trajectories yields the total spin current on 
the left side of the junction, l[{t) = {j[{pf,t)), where 

— (...) = /p^^^ (iV'sin(2V')(...)/4, and V = arccos(v/ ■ fi). 
Eqs. ([H])-® are our central results. They are valid for 
very general contact scattering matrices. 

In general spin-active junctions do not conserve the spin 
current. However, for the SFS point contacts described 
above, only the z component of spin current is nonzero 
and it is conserved across the contact, Is{t) = Is{t)z- 
Furthermore, the suppression of order parameter near 
the point contact, as well as spin accumulation in the 



and 



leads, are negligible. Thus, 7j^(e), lf{^), xf{t), 
ij^(e) can be evaluated in terms of their bulk values in 
the corresponding leads ,19]. For the results reported be- 
low we neglected the directional dependence of and d 
and assumed an inelastic scattering rate of 10~^A. The 
IV characteristics, when normalized by the temperature- 
dependent gap A(r), depend only weakly on the tem- 
perature, so we only present results for T — Q. 

Figure [5] shows the dc component of spin current at zero 
temperature for two sets of spin dependent transparen- 
cies, and three values of the spin mixing angle. Con- 
sider first the case of pure spin filtering, i.e. ?? = 0. 
For eV > 2A, the spin current is always less than the 
normal state value, GsV, where the spin conductance 
Gs = NfVf^A{Di; - D{). The SGS of the spin cur- 
rent shows anomalies at voltages, Vn — 2A/ne, which 
are more pronounced in higher transmission junctions. 
The structures associated with even and odd harmonics 
are different: the spin current is a local minimum near 
V2m+i, but a local maximum near V2m- The SGS for 
the spin current is in sharp contrast with that of the 
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FIG. 3: Left (right): The cosine (sine) components of the ac 
spin current (in units of C?sA/e) at T" = for = 0, tt/2, tt, 
= 0.5, and = 0.4. 



charge current |16[. The deficit in spin current at high 
voltages, as well as the even-odd effect in the SGS reflect 
the quasiparticle origin of the spin current. Multiple An- 
dreev reflection processes starting at subgap energies con- 
tribute significantly to the charge current, but much less 
so to the spin current. Consider a scattering event where 
an incoming spin-up electron-like quasiparticle of energy 
e < — A in lead 1 climbs up the MAR ladder to energy 
e„ > A after n—1 Andreev reflections [n ~ 2 A /eV) and 
leaks into the leads. This process transfers charge ne to 
lead 2, but only spin angular momentum /i/2, and only 
if n is odd. The rfi^ order MAR process has an onset 
voltage, Vn- Opening up new transmission channels at 
V2m+i leads to increase of the spin current, while the on- 
set of back scattering at V2m leads to decrease in spin 
current. For larger contact transparencies contributions 
from higher-order MAR processes become relevant, and 
more detailed structures in the IV curve develop at low 
voltages, as are evident in the right panel of Fig. [2l 

The inclusion of spin mixing leads to interface bound 
states at energies Acos(?9/2) [spin up] and — Acos(?9/2) 
[spin down] [l5|. Qualitatively, this structure of the lo- 
cal density of states (LDOS) affects the spin transport 
in two ways. First, the LDOS for continuum states near 
the gap edge is depleted. This leads to suppression of 
the spin current at high voltages; also the SGS at the on- 
set voltages for MAR processes is less pronounced. Sec- 
ond, scattering events with MAR trajectories which coin- 
cide with the bound states exhibit resonant transmission 
and dominate the spin current. Thus, in relatively low- 
transmission junctions the SGS is displaced from Vn to 
voltages near A(l -I- cos(i9/2))/n Ti*]. Since the reso- 
nances for spin up and spin down quasiparticles occur 
at different voltages, the dc spin current oscillates with 
voltage, even reversing direction in junctions with mod- 
erate to high transparency. For i9 = tt, the bound states 
are degenerate and located at the Fermi level, and con- 
sequently the SGS at Vn is recovered. 

Figure[3]shows the cosine and sine part of the first Fourier 
component of the ac spin current at zero temperature for 



intermediate transparencies: = 0.5, Di = 0.4. Note 
that the magnitudes of /f and If are comparable to that 
of the dc current. However, higher order Fourier com- 
ponents of the ac spin current are generally smaller in 
magnitude, but readily calculated from Eqs. ([S])-®- As 
in the case of the dc current, the structures in Ii{V) 
and IfiV) are related to onsets and resonances associ- 
ated with MAR processes. Resonant transmission leads 
to dramatic variations and sharp peaks in the ac spin 
current at voltages less than 2A; e.g. for = tt the sine- 
component dominates the response near eV/A = 0.5 and 
reverses sign over a narrow voltage range, eSV/ A « 0.1. 

In conclusion, we calculated the dynamical spin current 
in voltage-biased spin-active Josephson contacts. The de- 
pendence of the spin current on the voltage bias exhibits 
nonlinearities and resonant structures that are in striking 
contrast to the IV characteristics of the charge current, 
but which can be understood in terms of spin-filtering, 
spin-mixing and MAR. The nonlinear spin IV charac- 
teristics may be used to control the magnitude and the 
direction of the spin current via bias voltage, and may 
be of interest for superconducting spintronic devices. 

We thank Tomas Lofwander for discussions on ac Joseph- 
son effects. 
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